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ABSTRACT 

We analyze high signal-to-noise spectrophotonietric observations acquired simul- 
taneously with TWIN, a double-arm spectrograph, from 3400 to 10400 A of three 
star-forming regions in the Hii galaxy SDSS J165712. 754-321141. 4. We have measured 
four line temperatures: Te([Oiii]), Te([Siii]), Te([Oii]), and Te([Sii]), with high preci- 
sion, rms errors of order 2%, 5%, 6% and 6%, respectively, for the brightest region, and 
slightly worse for the other two. The temperature measurements allowed the direct 
derivation of ionic abundances of oxygen, sulphur, nitrogen, neon and argon. 

We have computed CLOUDY tailor-made models which reproduce the 0^+ mea- 
sured thermal and ionic structures within the errors in the three knots, with deviations 
of only 0.1 dex in the case of 0+ and S^+ ionic abundances. In the case of the electron 
temperature and the ionic abundances of S^/H"*", we find major discrepancies which 
could be consequence of the presence of colder diffuse gas. The star formation history 
derived using STARLIGHT shows a similar age distribution of the ionizing popula- 
tion among the three star-forming regions. This fact suggests a similar evolutionary 
history which is probably related to a process of interaction with a companion galaxy 
that triggered the star formation in the different regions almost at the same time. 
The hardness of the radiation field mapped through the use of the softness parame- 
ter rj is the same within the observational errors for all three regions, implying that 
the equivalent effective temperature of the radiation fields are very similar for all the 
studied regions of the galaxy, in spite of some small differences in the ionization state 
of different elements. 

Regarding the kinematics of the galaxy, the gas rotation curve shows a deviation 
from the circular motion probably due either to an interaction process or related to 
an expanding bubble or shell of ionized gas approaching us. A dynamical mass of 
2.5 X 10"'^° M0 is derived from the rotation curve. 

Key words: ISM: abundances - Hii regions - galaxies: abundances - galaxies: funda- 
mental parameters - galaxies: starburst - galaxies: stellar content. 



1 INTRODUCTION 
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Star formation is an ongoing process in the local universe, 
with observed rates of the order of lO"'^ Mq yr"'^ Mpc~^ 
(Madau et al. 1996). Most of the light and metals are pro- 
duced in the most massive among the newly formed stars. 
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The most extreme regions forming massive stars are often 
referred to as starbursts. In the local universe they account 
for about a quarter of all star formation (Heckman 1997), 
and this fraction may have been larger in the younger uni- 
verse. The origin of the term "starburst" (coined as "star- 
burst nuclei" by Weedman et al. 1981) dates back to the 
early observations of dust-obscured star-forming regions in 
the centres of nearby galaxies at the end of the seventies 
and beginning of the eighties, but the basic concept extends 
further back {e.g., Hodge 1969; Searle et al. 1973). 

The level of intensity of a starburst is highly variable. 
According to Terlevich (1997), in a starburst galaxy the en- 
ergy output of the starburst (Lsb) is much larger than the 
one coming from the rest of the galaxy (Lg), a galaxy with 
LsB ~ Lg is a galaxy with starbursts, and in a normal 
galaxy Lsb ^ Lg. This classification shows the variety of 
environments of the bursts. It is clear that the visibility of 
the burst depends not only on its intensity but also on its 
environment. Terlevich (1997) also proposed a division in 
phases of the starburst. The first one, the nebular phase, is 
characterized by the presence of strong emission lines from 
gas photoionized by young massive stars, with an age of less 
than 10 Myr. The early continuum phase goes from 10 to 100 
Myr, when some Balmer lines appear in absorption and oth- 
ers in emission. Finally, the late continuum phase, is when 
only some weak emission lines appear in the spectrum. The 
Hll galaxies are typical examples of the first phase. 

Hll galaxies are gas-rich dwarf galaxies experiencing a 
violent star formation period which dominates the optical 
spectrum of the host galaxy. They have one of the high- 
est intensity levels of star forming activity. In general, these 
galaxies have a central region which contains one or more 
star forming knots, with a diameter of several hundred par- 
sees with high surface brightness, and a low luminosity un- 
derlying galaxy (My ^ -17). The activity of the star for- 
mation episodes can not be sustained continuously for long 
periods of time, since the central region can not have enough 
gas to fuel these processes for longer than 10^ years and to 
match the gas content and metallicity with theoretical con- 
siderations (Thuan et al. 2004). 

Spectroscopically, Hll galaxies are essentially identical 
to the giant Hll regions found in nearby irregular and late- 
type galaxies. The correlation among structural parameters 
(H/3 luminosity, velocity dispersion, line widths) and be- 
tween these parameters and chemical composition (Terlevich 
& Melnick 1981) favours the interpretation of Hll galaxies as 
giant Hll regions in distant dwarf irregular galaxies similar 
to the ones found nearby (Melnick et al. 1985). 

Other important characteristic of Hll galaxies is their 
low metallicity (Z0/5O ^5 Z ^ Zq/S; Kunth & Sargent 
1983). The fact that Hll galaxies are metal-poor and very 
blue objects seems to suggest that they are young. Nev- 
ertheless, there is evidence which indicates the presence of 
populations older than the ones in the starburst. This is 
seen in the behaviour of the surface brightness profile which 
is exponential in the external zones, or in the colour index, 
which turns redder in V-R and V-I (Telles & Terlevich 1997). 
IZwl8 in particular, was considered as the best candidate 
for a truly young galaxy. Early studies of the stellar popula- 
tion of IZwl8 did not reveal any old population (Hunter & 
Thronson 1995). This contradicted some models which pre- 
dict that during a starburst, the heavy elements produced 



by the massive stars are ejected with high velocities into a 
hot phase, leaving the starburst region without immediate 
contribution to the enrichment of the interstellar medium 
(Tenorio-Tagle 1996). In this scenario, the metals observed 
now would have their origin in a previous star formation 
event, and an underlying old stellar population would be 
expected. In fact, Garnett et al. (1997) attributed the high 
carbon abundance that they found in HST spectroscopy of 
IZw 18 as evidence for the presence of an old stellar popula- 
tion. In agreement with this result, using HST archive data 
(Aloisi et al. 1999) showed that stars older than 1 Gyr must 
be present in IZw 18. Moreover, studies of the resolved stellar 
population in the near infrared with NICMOS (Osthn 2000) 
found also that while the near infrared colour-magnitude 
diagram was dominated by stars 10-20 Myr old, the pres- 
ence of numerous AGB stars require an age of at least 10* 
years. Legrand et al. (2000) modelled the relative abundance 
of metals in IZw 18 and concluded that, in addition to the 
present burst of star formation, a low star formation rate ex- 
tended over a long period of time was necessary to account 
for the observed values. 



In recent years, with the development of the Integral 
Field Unit (IFU) instruments to perform 3D spectroscopy, 
works that require a spatial coverage to study extended 
galactic or extra galactic star-forming regions have been 
mainly focused on the use of this technique (Relano et al. 
2010; Cairos et al. 2010; Monreal-Ibero et al. 2010; Rosales- 
Ortega et al. 2010; Perez-Gallego et al. 2010; Garci'a-Benito 
et al. 2010; Sanchez et al. 2010; Perez-Montero et al. 2011, 
see for example). However, medium or high dispersion slit 
spectroscopy are a better option for spectrophotometry, or 
when the object is very compact, or even extended but 
with few star-forming knots. This is also the case when 
good spatial and spectral resolution and simultaneous wide 
spectral coverage are required (Diaz et al. 2007; Gumming 
et al. 2008; Firpo et al. 2010; Hagele et al. 2006, 2007, 
2008; Hagele 2008; Hagele et al. 2009, 2010; Perez-Montero 
et al. 2009; Lopez-Sanchez & Esteban 2009, 2010a,b; Lopez- 
Sanchez 2010; Firpo et al. 2011, see for example). 



In this paper we present simultaneous blue and red 
long-slit observations obtained with the double arm TWIN 
spectrograph at the 3.5m telescope of Calar Alto of the 
three brightest star-forming knots of the Hll galaxy SDSS 
J165712. 75+321141.4. This is part of a project to obtain 
a top quality spectrophotometric data base to determine 
ionized gas parameters which are indispensable to critically 
test photoionization models and to explore discrepancies be- 
tween models and observations. In Section 2 we show the 
details of the observations and data reduction. Section 3 
presents the derived physical characteristics of the regions, 
including the electron temperature for four different species. 
Section 4 is devoted to the discussion of these results, and fi- 
nally the summary and conclusions are presented in Section 
5. 
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Table 1. Right ascension, declination, redshift and SDSS photometric magnitudes of the observed knots obtained using the SDSS explore 
tools". 



Object ID 
(spSpcc SDSS) 


hereafter ID 


Knot 




RA 




Dec 


redshift 


u 


g 


r 


i 


z 


SDSS J165712. 75+321141.4 


SDSS J1657 


A 


le*" 


57'" 12! 75 


32° 


11'41'.'42 


0.038 


17.63 


17.01 


17.25 


17.14 


17.16 


(spSpec-52791-1176-591) 




B 


le'' 


57'" 12?26 


32° 


11' 43'.' 20 




20.34 


19.76 


20.65 


20.72 


20.27 






C 


16'' 


57"^ 13?58 


32° 


11' 40'.'09 




19.28 


18.75 


19.52 


19.40 


19.17 



"http://cas. sdss.org/ astro/en/ tools / explore /ob j . asp 




2 OBSERVATIONS AND DATA REDUCTION 

2.1 Object selection 

Using the implementation of the SDSS database in the 
INAOE Virtual Observatory superserver^, we selected the 
brightest nearby narrow emission line galaxies with very 
strong lines and large equivalent widths of Hq from the 
whole SDSS data release available at the time of planning 
the observations. These preliminary lists were then pro- 
cessed using BPT (Baldwin, Phillips & Terlevich 1981) di- 
agnostic diagrams to remove AGN-like objects. The final 
list consisted of about 10500 bonafide bright Hll galaxies. 
They show spectral properties indicating a wide range of 
gaseous abundances and ages of the underlying stellar pop- 
ulations (Lopez 2005). From this list, the final set was se- 
lected by further restricting the sample to the largest Ha 
flux and highest signal-to-noise ratio objects (for a complete 
description of the selection criteria see Hagele et al. 2006, 
hereafter Paper I). Of the selected sample, we chose SDSS 
J165712. 75-1-321141.4 to be observed at the allocated time. 
For simplicity, we will call the galaxy SDSS J1657 in what 
remains of the paper. 

Some general characteristics of the knots of SDSS J1657 
collected from the SDSS web page are listed in Table 1. 

2.2 Observations 

Blue and red spectra were obtained simultaneously using 
the double beam Cassegrain Twin Spectrograph (TWIN) 
mounted on the 3.5m telescope of the Calar Alto Observa- 
tory at the Centro Astronomico Hispano Aleman (CAHA), 
Spain. These observations were part of a four night observing 
run in 2006 June and they were acquired under excellent see- 
ing and photometric conditions (for details see; Hagele et al. 
2008, hereinafter Paper II). The blue arm covers the wave- 
length range 3400-5700 A (centred at Ac = 4550 A), giving 
a spectral dispersion of 1.09Apixel~^ (R~4170). On the 
red arm, the spectral range covers from 5800 to 10400 A 
(Ac = 8100 A) with a spectral dispersion of 2.42Apixel~^ 
(R~3350). The slit width was ~1.2arcsec. The pixel size 
for this set-up configuration is 0.56 arcsec for both spectral 
ranges. The target was observed at paralactic angle to avoid 
effects of differential refraction in the UV. As it can be seen 
in Fig. 1, the three main knots of SDSS J1657 are almost per- 
fectly aligned along the paralactic angle. The instrumental 
configuration, summarized in Table 2, covers the whole spec- 
trum from 3400 to 10400 A (with a gap between 5700 and 
5800 A) providing a moderate spectral resolution. This spec- 
tral coverage guarantees the simultaneous detection of the 

^ http:/ /astro. inaoep.mx/en/observatories/virtual/ 



Figure 1. False colour image of SDSS J1657 with the slit po- 
sition and the adopted knot names superimposed. This image 
was obtained using the SDSS explore tools. Circles and squares 
represent the photometric and spectroscopic SDSS targets, re- 
spectively. The scale is 782 pc arcsec"^ , at the adopted distance 
for SDSS J1657. [See the electronic edition of the Journal for a 
colour version of this figure.] 



Table 2. CAHA instrumental configuration. 





Spectral range 


Disp. 


^FWHM 


Spatial res. 




(A) 


(Apx-i) 


("px-l) 


Blue 


3400-5700 


1.09 


1420 


0.56 


Red 


5800-10400 


2.42 


1160 


0.56 



"RpWHM = -^/^^FWHM 

nebular fines from [On] AA 3727,29 to [Siii] AA 9069,9532 A 
at both ends of the spectrum, in the very same region of the 
galaxy, with a good signal-to- noise ratio (S/N) that allows 
the measurement of the weak auroral lines. Typical values of 
S/N are ~60 for [Oin]A4363 and ~20 for [Sii]A4068 (see 
Table 4 of Paper II). 



2.3 Data reduction 

Several bias and sky flat field frames were taken at the be- 
ginning and at the end of the night in both arms. In addi- 
tion, two lamp flat flelds and one He-Ar calibration lamp 
exposures were performed at each telescope position. The 
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-10 

distance (arcsecs) 



Figure 2. Spatial profile of the light distribution along the 
slit for the observed Ho emission. The profiles correspond to 
line+continuum (dashed line), continuum (dashed-dottcd line) 
and the difference between them (solid line) , representing the pure 
emission from Ha. 



images were processed and analysed with IRAF'^ routines 
in the usual manner. This procedure includes the removal 
of cosmic rays, bias subtraction, division by a normalised 
flat field, and wavelength calibration. To finish, the spectra 
are corrected for atmospheric extinction and flux-calibrated. 
Four standard star observations were performed each night 
at the same time for both arms, allowing a good spectropho- 
tometric calibration with an estimated rms error of about 
3%. Further details concerning each of these steps can be 
found in Paper II. 

Fig. 2 shows the spatial distribution of the Ha flux and 
the continuum along the slit for SDSS J1657. The emission 
line profiles have been generated by collapsing 11 pixels of 
the spectra in the direction of the resolution at the central 
position of the emission lines in the rest frame, A 6563 A, 
and are plotted as a dashed line. Continuum profiles were 
generated by collapsing 1 1 resolution pixels centred at 30 A 
to the red for each region and are plotted as a dashed-dotted 
line. The difference between the two, shown as a solid line, 
corresponds to the pure emission. From Fig. 2 it is clear that 
the continuum emission is not very strong, specially in the 
weaker knots. The three regions are labeled in the figure. 
There is a weak pure emission knot located between Knot 
A and C, which does not have enough S/N to derive the 
physical conditions of the gas. 



^ IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatories, which is oper- 
ated by the Association of Universities for Research in Astronomy, 
Inc. (AURA) under cooperative agreement with the National Sci- 
ence Foundation (NSF). 



3 RESULTS 

3.1 Line intensities and reddening correction 

The spectra of the three knots of SDSS J1657 (labelled from 
A to C) with some of the relevant identified emission lines 
are shown in Fig. 3. The spectrum of each observed knot is 
split into two panels. Knot A corresponds to the one anal- 
ysed in Paper II. 

The emission line fluxes were measured using the splot 
task in iraf following the procedure described in Paper I. 
Following Perez-Montero & Diaz (2003), the statistical er- 
rors associated with the observed emission fluxes have been 
calculated using the expression 

m = (JciV'''^[l4-£W/(iVA)]'''^ 

where ai is the error in the observed line flux, ctc represents 
the standard deviation in a box near the measured emission 
line and stands for the error in the continuum placement, N 
is the number of pixels used in the measurement of the line 
flux, EW is the line equivalent width, and A is the wave- 
length dispersion in A per pixel (Gonzalez-Delgado et al. 
1994). There are several emission lines affected by cosmetic 
faults or charge transfer in the CCD, internal reflections in 
the spectrograph, telluric emission lines or atmospheric ab- 
sorption lines. These cause the errors to increase, and, in 
some cases, they are impossible to quantify, in which case 
they were ignored and excluded from any subsequent anal- 
ysis. 

Some observed lines (e.g., [Cllll] AA 5517,5537, several 
carbon recombination lines, Balmer or Paschen lines) were 
impossible to measure due to low signal to noise. This is also 
the case for the Balmer and Paschen jump, that could not 
be measured because of the difficulty to fit the continuum at 
both sides of the discontinuity with an acceptable precision. 

An underlying stellar population is easily appreciable by 
the presence of absorption features that depress the Balmer 
and Paschen emission lines. A pseudo-continuum has been 
defined at the base of the hydrogen emission lines to mea- 
sure the line intensities and minimize the errors introduced 
by the underlying population (see Paper I). The presence of 
the wings of the absorption lines imply that, even though 
we have used a pseudo-continuum, there is still an absorbed 
fraction of the emitted flux that we are not able to measure 
accurately (see discussion in Di'az 1988). This fraction is not 
the same for all lines, nor are the ratios between the ab- 
sorbed fractions and the emission. In Paper I we estimated 
that the difference between the measurements obtained us- 
ing the defined pseudo-continuum or a multi- Gaussian fit 
to the absorption and emission components is, for all the 
Balmer lines, within the errors. This is also the case for 
the objects studied here. At any rate, for the Balmer and 
Paschen emission lines we have doubled the derived error, 
0";, as a conservative approach to include the uncertainties 
introduced by the presence of the underlying stellar popu- 
lation. 

The absorption features of the underlying stellar pop- 
ulation may also affect the helium emission lines to some 
extent. However, these absorption lines are narrower than 
those of hydrogen (see, for example, Gonzalez-Delgado et al. 
2005). Therefore it is difficult to set adequate pseudo- 
continua at both sides of the lines to measure their fiuxes. 
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Figure 3. Blue and red spectra of Knots A, B and C of SDSS J1657 (upper, middle and lower panel, respectively) in the rest frame. 
The flux scales are the same in both spectral ranges. 



We also applied the STARLIGHT code^ (Cid Fernan- 
des et al. 2005) to each region to separate the emission 
spectra from the underlying stellar absorptions, but for the 
strongest emission lines the difference between the measure- 
ments made after the subtraction of the STARLIGHT fit 
and the ones made using the pseudo-continuum is well be- 
low the observational errors. For a detailed discussion on 
the differences in the emission line measurements see Perez- 
Montero et al. (2010). 



^ The STARLIGHT project is supported by the Brazilian agen- 
cies CNPq, CAPES and FAPESP and by the France-Brazil 
CAPES/Cofecub program. 



The reddening coefficients c(H/3) were calculated 
from the measured Balmer decrements, F{X) / F{H P). We 
adopted the galactic extinction law of Miller & Mathews 
(1972) with 7?v=3.2. A least square fit of the measured 
decrements to the theoretical ones, {F{X) / F{H l3))o, com- 
puted based on the data by Storey & Hummer (1995), was 
performed that provides the value of c(H/3). The theoretical 
Balmer decrements depend on electron density and temper- 
ature. We used an iterative method to estimate them, taking 
as starting values those derived from the measured [Sll] AA 
6717,6731 A and [Oiii] AA 4363, 4959, 5007 A line fluxes. Due 
to the large error introduced by the presence of the under- 



6 G. F. Hdgele et al. 



Table 3. Relative reddening corrected line intensities [-F(i?/3)=/(_ff/3)=10000] for the three star-forming knots. 



A (A) 






Knot A 






Knot B 






Knot C 




f(A) 


EW 


I{X) 


Error 


EW 


/(A) 


Error 


EW 


/(A) 


Error 






(A) 




(%) 


(A) 




(%) 


(A) 




(%) 


3727 [Oii]° 


0.271 


120.1 


18832 ± 230 


1.2 


107.8 


13266 ±420 


3.2 


112.2 


14809 ±282 


1.9 


3750 H12 


0.266 


1.9 


232 ± 47 


20.2 








_ 


_ 




3770 Hll 


0.261 


2.3 


293 ± 40 


13.8 


9.7 


617 ±146 


23.7 











3798 HIO 


0.254 


4.1 


500 ± 68 


13.5 








4.9 


465 ± 141 


30.4 


3835 H9 


0.246 


7.1 


780 ± 93 


11.9 


8.9 


783 ± 180 


23.0 


7.5 


761 ± 140 


18.4 


3868 [Neiii] 


0.238 


23.1 


3262 ±132 


4.0 


30.7 


3663 ±132 


3.6 


28.6 


3533 ±118 


3.3 


3889 Hei+H8 


0.233 


14.0 


1826 ± 95 


5.2 


22.3 


1954 ± 208 


10.7 


26.2 


2192 ±267 


12.2 


3968 [Neiii]+H7 


0.216 


22.4 


2456 ±121 


4.9 


30.1 


2705 ±221 


8.2 


29.9 


2699 ± 237 


8.8 


4026 [Nll]+Hel 


0.203 


1.1 


155 ± 16 


10.4 


4.3 


327 ±105 


32.1 








4068 [Sii] 


0.195 


1.4 


198 ± 15 


7.4 


1.4 


153 ± 36 


23.7 


1.2 


130 ± 35 


26.9 


4102 H5 


0.188 


20.9 


2432 ± 65 


2.7 


30.2 


2667 ±163 


6.1 


29.2 


2619 ±173 


6.6 


4340 H7 


0.142 


43.1 


4417 ± 97 


2.2 


76.5 


4813 ±219 


4.5 


68.7 


4712 ± 140 


3.0 


4363 [Oiii] 


0.138 


4.5 


524 ± 24 


4.6 


8.3 


846 ± 67 


7.9 


8.9 


838 ± 59 


7.0 


4471 Hei 


0.106 


4.2 


443 ± 33 


7.4 


4.7 


405 ± 44 


11.0 


5.3 


445 ± 41 


9.3 


4658 [Fein] 


0.053 


1.0 


107 ± 16 


14.9 














4686 Hen 


0.045 


1.2 


126 ± 14 


11.0 


4.3 


313 ± 68 


21.9 


2.1 


170 ± 24 


14.2 


4861 H/3 


0.000 


117.8 


10000 ± 79 


0.8 


153.3 


10000 ± 178 


1.8 


167.2 


10000 ±128 


1.3 


4921 Hei 


-0.014 


0.8 


75 ± 14 


18.6 














4959 [Oni] 


-0.024 


152.5 


14333 ± 127 


0.9 


218.5 


16118 ±129 


0.8 


199.2 


14940 ±114 


0.8 


4986 [Fem]'' 


-0.030 


1.4 


135 ± 28 


20.5 


2.4 


163 ± 57 


35.1 


3.2 


215 ± 39 


18.3 


5007 [Oni] 


-0.035 


455.1 


43082 ± 240 


0.6 


705.5 


48653 ± 256 


0.5 


613.2 


44727 ± 129 


0.3 


5015 Hei 


-0.037 


2.4 


222 ± 23 


10.1 


3.2 


205 ± 35 


17.3 


3.2 


220 ± 28 


12.6 


5199 [Ni] 


-0.078 


2.0 


157 ± 26 


16.4 














5876 Hei 


-0.209 


18.9 


1116 ± 44 


3.9 


23.4 


991 ± 29 


2.9 


29.0 


1149 ± 47 


4.1 


6300 [Oi] 


-0.276 


8.1 


438 ± 16 


3.7 


9.6 


393 ± 26 


6.7 


8.6 


366 ± 14 


3.9 


6312 [Sin] 


-0.278 


3.7 


201 ± 9 


4.3 


4.3 


175 ± 10 


5.8 


3.5 


148 ± 10 


6.4 


6364 [Oi] 


-0.285 


2.8 


152 ± 18 


11.8 


3.6 


141 ± 20 


14.4 


3.0 


124 ± 12 


10.0 


6548 [Nil] 


-0.311 


9.4 


464 ± 23 


5.0 


5.5 


216 ± 23 


10.7 


6.7 


266 ± 14 


5.4 


6563 Ho 


-0.313 


571.3 


27772 ± 153 


0.5 


772.5 


28159 ±105 


0.4 


730.5 


27919 ±133 


0.5 


6584 [Nil] 


-0.316 


28.8 


1428 ± 47 


3.3 


16.2 


632 ± 50 


8.0 


17.0 


680 ± 29 


4.3 


6678 Hei 


-0.329 


6.7 


315 ± 18 


5.7 


7.3 


272 ± 11 


4.1 


7.9 


296 ± 20 


6.9 


6717 [Sii] 


-0.334 


47.4 


2207 ± 57 


2.6 


39.4 


1489 ± 35 


2.3 


45.1 


1699 ± 56 


3.3 


6731 [Sii] 


-0.336 


32.2 


1598 ± 43 


2.7 


28.1 


1060 ± 25 


2.3 


30.7 


1154 ± 40 


3.4 


7065 Hei 


-0.377 


5.6 


235 ± 10 


4.4 


7.5 


226 ± 12 


5.1 


8.8 


279 ± 26 


9.3 


7136 [Ann] 


-0.385 


16.4 


717 ± 26 


3.6 


18.4 


584 ± 19 


3.3 


17.6 


581 ± 17 


2.9 


7281 Hei= 


-0.402 


0.9 


41 ± 7 


18.2 


3.0 


90 ± 14 


15.9 


1.4 


46 ± 13 


28.6 


7319 [On]'' 


-0.406 


12.3 


302 ± 17 


5.6 


4.8 


165 ± 15 


9.2 


6.1 


196 ± 18 


9.2 


7330 [011]"= 


-0.407 


8.8 


211 ± 14 


6.5 


7.2 


251 ± 21 


8.4 


7.9 


254 ± 26 


10.1 


7751 [Arm] 


-0.451 


4.6 


177 ± 22 


12.3 


5.5 


163 ± 13 


8.0 


5.8 


170 ± 17 


9.9 


8665 P13 


-0.531 


7.7 


144 ± 53 


37.1 








8.9 


133 ± 38 


28.3 


8751 P12 


-0.537 


4.2 


101 ± 29 


28.2 














8865 Pll 


-0.546 


8.5 


211 ± 34 


16.3 


17.9 


277 ± 73 


26.2 


11.3 


228 ± 46 


20.2 


9014 PIO 


-0.557 


15.4 


167 ± 36 


21.5 














9069 [Sin] 


-0.561 


59.2 


1400 ± 99 


7.1 


68.2 


984 ±122 


12.4 


81.9 


1577 ±134 


8.5 


9229 P9 


-0.572 


16.9 


263 ± 47 


18.0 


26.1 


364 ±129 


35.6 








9532 [Siii] 


-0.592 


157.3 


3674 ± 257 


7.0 


85.1 


2700 ± 152 


5.6 


238.7 


2915 ±162 


5.6 


I(H/3)(erg sec~-^ 


cm~"^) 




6.3 X 10"^'' 






2.6 X lO"^'' 






3.4 X lO"^'' 





c(H/3) 0.05 ±0.01 0.15 ±0.02 0.13 ±0.02 

" [Oil] AA 3726 + 3729; " [Fciii] AA 4986 + 4987; " possibly blend with an unknown line; " [On] AA 7318 + 7320; [On] AA 7330 + 7331. 



lying stellar population, only the strongest Balmer emission 
lines (Hof, H/3, H7 and H(5) were used. 

For the easiness of comparison, we have included in the 
following sections the results presented in Paper II for knot 
A. Table 3 lists the reddening corrected emission lines for 
each knot, together with the reddening constant and its er- 
ror taken as the uncertainties of the least square fit and the 
reddening corrected II/3 intensity. Column 1 shows the wave- 
length and the name of the measured lines. The adopted 
reddening curve, /(A), normalized to H/3, is given in col- 
umn 2. The errors in the emission lines were obtained by 



propagating in quadrature the observational errors in the 
emission line fluxes and the reddening constant uncertain- 
ties. We have not taken into account errors in the theoretical 
intensities since they are much lower than the observational 
ones. 



3.2 Physical conditions of the gas 

The physical conditions of the ionized gas, including electron 
temperatures (Te) and electron density (Ne ~n([Sll])), have 
been derived from the emission line data using the same 
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Table 4. Electron densities and temperatures. Densities in cm ^ 
and temperatures in 10* K. 

n([Sii]) Te([Oni]) Te([On]) Te([Sin]) T,([Sn]) 

1.23±0.02 1.33±0.07 1.45±0.08 0.88±0.05 

1.43±0.05 1.52±0.12 1.64±0.11 1.00±0.17 

1.48±0.05 1.50±0.13 1.29±0.07 0.83±0.15 



Knot A 


30 


Knot B 


10 


Knot C 


10 



Table 5. Ionic and total helium abundance. 



Knot A 



Knot B 



Knot C 



He+/H+(A4471) 

He+/H+(A5876) 

He+/H+(A6678) 

He+/H+(A7065) 

He+/H+(Adop.) 

He2+/H+(A4686) 

(He/H) 



0.093±0.007 
0.085±0.003 
0.086±0.005 
0.093±0.005 
0.087±0.005 
O.OOllitO.OOOl 
0.088±0.008 



0.085±0.009 
0.079±0.002 
0.077±0.003 
0.086±0.006 
0.080±0.005 
0.0028±0.0006 
0.080±0.008 



0.094±0.009 
0.093±0.003 
0.084±0.005 
0.105±0.010 
0.092±0.009 
0.0015±0.0002 
0.092±0.009 



procedures as in Paper II, based on the five-level statistical 
equilibrium atom approximation in the task temden, of the 
software package IRAF (de Robertis et al. 1987; Shaw & 
Dufour 1995). As usual, we have taken as sources of error 
the uncertainties associated with the measurement of the 
emission-line fluxes and the reddening correction, and we 
have propagated them through our calculations. 

For all three knots we have derived the elec- 
tron temperatures of [Oil], [Olll], [Sll] and [Slll]. The 
[Oil] AA 7319,7330 A lines have a contribution by direct re- 
combination which increases with temperature. Using the 
calculated [Olll] electron temperatures, we have estimated 
these contributions to be less than 4 % in all cases and there- 
fore we have not corrected for this effect. The expression for 
the correction of direct recombination, however, is valid only 
in the range of temperatures between 5000 and 10000 K. 
The temperatures found are slightly over that range. At any 
rate, the relative contribution of recombination to coUisional 
intensities decreases rapidly with increasing temperature, 
therefore for the high Te values found in our objects this 
contribution is expected to be small. 

The derived electron densities and temperatures for the 
three star-forming regions are given in Table 4 along with 
their corresponding errors. 



3.3 Chemical abundance derivation 

We have derived the ionic chemical abundances of the dif- 
ferent species using the strongest available emission lines 
detected in the analyzed spectra and the task ionic of the 
STSDAS package in IRAF, as described in Paper II. 

The total abundances have been calculated by taking 
into account, when required, the unseen ionization stages 
of each element, using the appropriate ionization correction 
factor (ICF) for each species, X/H = ICF(X+') X+7H+ as 
detailed in what follows. 



3. 3. 1 Helium 

We have used the well detected Hei AA 4471, 5876, 6678 and 
7065 A lines, to calculate the abundances of once ionized he- 
lium. For the three knots also the Hen A 4686 A line was mea- 
sured allowing the calculation of twice ionized He. The He 



lines arise mainly from pure recombination, although they 
could have some contribution from collisional excitation and 
be affected by self-absorption (see Olive & Skillman 2001, 
2004, for a complete treatment of these effects). We have 
taken the electron temperature of [Olll] as representative of 
the zone where the He emission arises since at any rate ra- 
tios of recombination lines are weakly sensitive to electron 
temperature. We have used the equations given by Olive & 
Skillman to derive the He'''/H^ value, using the theoretical 
emissivities scaled to H/? from Benjamin et al. (1999) and the 
expressions for the collisional correction factors from King- 
don & Ferland (1995). To calculate the abundance of twice 
ionized helium we have used equation (9) from Kunth & 
Sargent (1983). We have not made any corrections for flu- 
orescence since three of the used helium lines have a small 
dependence with optical depth effects but the observed ob- 
jects have low densities. We have not corrected either for 
the presence of an underlying stellar population. A sum- 
mary of the equations used to calculate these ionic abun- 
dances is given in Appendix B of Garcia- Benito (2009). The 
total abundance of He has been found by adding directly 
the two ionic abundances, He/H = (He+-hHe2+)/H+. The 
results obtained for each line and the total He abundances, 
along with their corresponding errors are presented in Table 
5. Also given in the table is the adopted value for He"'"/H"'" 
as the average, weighted by the errors, of the abundances 
derived from each Hei emission line . 



3.3.2 Ionic and total chemical abundances from forbidden 
lines 

The oxygen ionic abundance ratios, O^/H""" and O^^/H"*", 
have been derived from the [On] AA 3727,29 A and 
[Olll]AA4959, 5007 A lines respectively using for each ion 
its corresponding temperature. At the temperatures derived 
here, most of the oxygen is in the form of O^ and O'^"'', there- 
fore the approximation 0/H= (0^-|-0'^^)/H"'" is a valid one. 

S'''/H^ and S^^/H+, abundances have been derived 
using Te([Sll]) and Te([Slll]) and the fluxes of the 
[Sii] AA 6717,6731 A and the near-IR [Siii] AA 9069, 9532 A 
emission lines, respectively. Unlike oxygen, a relatively im- 
portant contribution from S'^"'' may be expected for sul- 
phur depending on the nebular excitation. The total sulphur 
abundance has been calculated using an ICF for S'*"-fS^''' 
according to the formula by Barker (1980), which is based 
on Stasihska (1978) photo-ionization models, with q = 2.5, 
which gives the best fit to the scarce observational data on 
S^^ abundances (Perez-Montero et al. 2006). Taking this 
ICF as a function of the ratio 0^'''/0 instead of O^/O re- 
duces the propagated error for this quantity. 

The ionic abundance of nitrogen, N^/H"*" has been de- 
rived from the intensities of the [Nil] AA 6548,6584 A lines as- 
suming Te([Nll])?aTe([Oll]). The N/O abundance ratio has 
been derived under the assumption that N/O = N"'"/©"'" and 
N/H was calculated as log(N/H) = log(N/0)-hlog(0/H). 

Neon is only visible in the spectra via the [Nelll] emis- 
sion line at A3868A, so Ne^"*" has been derived using this 
line. For this ion we have taken the electron tempera- 
ture of [Olll], as representative of the high excitation zone 
(Te([Nelll])«Te([Olll]); Peimbert & Costero 1969). Classi- 
cally, the total abundance of neon has been calculated as- 
suming that Ne/0«Ne^+/0^+. Izotov et al. (2004) point 
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Table 6. Ionic chemical abundances derived from forbidden emis- 
sion lines, ICFs and total chemical abundances for elements heav- 
ier than helium. 





Knot 


A 


Knot 


B 


Knot 


C 


12 -1- log(0+/H+) 
12 + log(02+/H+) 


7.37 ± 


0.09 


7.03 ± 


0.12 


7.10 ± 


0.12 


7.87 ± 


0.02 


7.74 ± 


0.04 


7.67 ± 


0.04 


12 + log(0/H) 


7.99 ± 


0.04 


7.82 ± 


0.05 


7.78 ± 


0.06 


ICF(0+-H02+)* 


1.07 


1.06 


1.06 


12 + log(0/H)* 


8.02 ± 


0.04 


7.84 ± 


0.05 


7.80 ± 


0.06 


12 + log(S+/H+) 


6.07± 


0.08 


5.76 ± 


0.19 


6.02 ± 


0.24 


12 + log(S2+/H+) 


6.00 ± 


0.06 


5.79 ± 


0.07 


6.02 ± 


0.07 


ICF(S+ -1- S2+) 


1.32 ± 


0.05 


1.51 ± 


0.10 


1.38 ± 


0.09 


12 + log(S/H) 


6.46 ± 


0.07 


6.25 ± 


0.13 


6.46 ± 


0.16 


log(S/0) 


-1.53 ± 


0.08 


-1.57± 


0.14 


-1.32± 


0.17 


ICF(S+ -1- S2+)* 


1.07 


1.18 


1.14 


12 -f log(S/H)* 


6.37 ± 


0.07 


6.15 ± 


0.13 


6.18 ± 


0.15 


log(S/0)* 


-1.65 ± 


0.08 


-1.69± 


0.14 


-1.42 ± 


0.16 


12 -1- log(N+/H+) 


6.15± 


0.06 


5.68 ± 


0.10 


5.74 ± 


0.09 


12 + log(N/H) 


6.76 ± 


0.28 


6.47 ± 


0.39 


6.42 ± 


0.39 


log(N/0) 


-1.23 ± 


0.11 


-1.35 ± 


0.16 


-1.36± 


0.15 


ICF(N+)* 


3.39 


4.90 


4.27 


12 + log(N/H)* 


6.68 ± 


0.06 


6.37 ± 


0.10 


6.37 ± 


0.09 


log(N/0)* 


-1.34 ± 


0.07 


-1.47± 


0.11 


-1.43± 


0.11 


12 -1- log(Ne^+/H+) 


7.22 ± 


0.04 


7.06 ± 


0.06 


7.00 ± 


0.05 


ICF(Ne) 


1.09 ± 


0.00 


1.08 ± 


0.00 


1.08 ± 


0.01 


12 + log(Ne/H) 


7.25 ± 


0.04 


7.09 ± 


0.06 


7.04 ± 


0.05 


log(Ne/0) 


-0.74 ± 


0.06 


-0.73 ± 


0.08 


-0.74 ± 


0.08 


ICF(Ne2+)* 


1.26 


1.12 


1.17 


12 4- log(Ne/H)* 


7.32 ± 


0.04 


7.11 ± 


0.06 


7.17± 


0.05 


log(Ne/0)* 


-0.70 ± 


0.05 


-0.73 ± 


0.08 


-0.73 ± 


0.08 


12 -1- log(Ar^+/H+) 


5.49 ± 


0.06 


5.31 ± 


0.06 


5.49 ± 


0.06 


ICF(Ar2+) 


1.13± 


0.02 


1.23 ± 


0.06 


1.16 ± 


0.04 


12 + log(Ar/H) 


5.54± 


0.06 


5.39 ± 


0.06 


5.55 ± 


0.06 


log(Ar/0) 


-2.45 ± 


0.07 


-2.43 ± 


0.08 


-2.22 ± 


0.08 



ICF(Ar2+)* 1.15 1.17 1.17 

12 -I- log(Ar/H)* 5.55 ± 0.06 5.38 ± 0.06 5.56 ± 0.06 

log(Ar/0)* -2.47 ill 0.07 -2.46ifc 0.08 -2.24ili 0.08 

* ICFs and total abundances from photoionization models (see 
text). 



out that this assumption can lead to an overestimate of 
Ne/H in objects with low excitation, where the charge trans- 
fer between O^^ and H" becomes important. Thus, we have 
used the expression of this ICF given by (Perez-Montero 
et al. 2007). It is interesting to note, however, that given 
the high excitation of the observed objects there is no sig- 
nificant difference between the total neon abundance derived 
using this ICF and those estimated using the classic approx- 
imation. 

The only available emission lines of argon in the optical 
spectra of ionized regions correspond to Ar^+ and Ar^^ . The 
abundance of Ar^+ has been calculated from the measured 
[Arm] A 7136 A line assuming that Te([Arlll]) Te([Slll]) 
(Garnett 1992). [Ariv] was not detected in the spectra. 
The total abundance of Ar was hence calculated using the 
ICF(Ar'^^) derived from photo- ionization models by Perez- 
Montero et al. (2007). 

The ionic abundances with respect to ionized hydrogen 
of the elements heavier than helium, ICFs, total abundances 
and their corresponding errors are given in Table 6. 



3.4 Photoionization models of the observed 
regions 

Detailed tailor-made photoionization models were produced 
in order to ascertain the main properties of the ionizing stel- 
lar population and the ionized gas. The methodology is de- 
scribed in Perez-Montero et al. (2010) who study the bright- 
est knots of the Hll galaxies described in Papers I and II, in- 
cluding knot A in SDSS J1657. Here we describe the models 
for knots B and C, and compare them with the observations 
and with the results obtained for knot A in Perez-Montero 
et al. (2010). 

We have resorted to the photoionization code CLOUDY 
V. 06.02c (Ferland et al. 1998), using the equivalent width 
of II;S, after removing the underlying stellar population (i.e. 
the population younger than 10 Myr), the Ha luminosity 
and the intensities of [On] 3727 A, [Oiii] 4363 and 5007 A, 
[Sll] 6717 and 6731 A, and [Slll] 9069 and 9532 A relative to 
H/3. 

We have used for the photoionization models the same 
Starburst 99 stellar libraries as in the model fitting of the 
stellar population, described in §4.4, with the metallicity 
closest to the value measured in the gas-phase, Z = 0.004 
(= 1/5 Zq). We assummed a constant star formation his- 
tory which, according to Perez-Montero et al., gives the 
best agreement for the number of ionizing photons and the 
EW(H/3) corrected for underlying stellar population and 
dust absorption effects. A thick shell geometry and a con- 
stant density of 100 particles per cm'^ have been set as input 
conditions in all the models. To fit the observed properties, 
the distance to the ionizing source, the filling factor, the 
dust-to-gas ratio and the age of the stellar cluster were left 
as free parameters. 

One of the most important parameters in the correct 
modelling of ionized gas nebulae is the dust absorption fac- 
tor, fd, which gives the ratio between the number of ionizing 
photons emitted by the stellar cluster and the number of 
ionizing photons absorbed by the gas (Perez-Montero et al. 
2010). This factor must be taken into account in deriving 
properties of the cluster from hydrogen Balmer recombina- 
tion lines. It has been obtained in the best model, after an 
iterative method to fit the observed relative emission-line 
intensities and the corrected EW(H/3) and L(Ha). In Fig. 4 
we show the ratio between the intensities of the most repre- 
sentative observed and modelled emission lines for the three 
knots. Data for knot A have been taken from Perez-Montero 
et al. (2010). As we can see, the agreement results excellent 
for all involved [Olll] lines, with a deviation smaller than 
5% in all three knots. In the case of the [On] lines and 
[Sin] 9069 A it is better than 10%. The fitting of [Slll] at 
6312 A is a bit worse, with a 20% of disagreement in knots 
A and C, and 30% in B. The largest discrepancy is found for 
the [Sll] lines, from 30% of disagreement up to 65% in the 
case of 6717,6731 A in knot B. In Table 7, we compare the 
observed and modelled EW(H/?), corrected for the contribu- 
tion of the underlying stellar population. We also give the 
number of ionizing photons and other properties predicted 
by the individual models, such as the age of the ionizing 
cluster, filling factor, ionization parameter, dust-to-gas ra- 
tio and visual extinction. Regarding knot A, all quantities 
have been extracted from the model in Perez-Montero et al. 
(2010). As we can see, the agreement between observed and 
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Figure 4. Ratio between observed and modelled intensities of the most representative emission lines for each one of the star-formig 
knots. 



Table 7. Observed and model-predicted properties of the three 
studied regions 







Knot A 


Knot B 


Knot C 


Age (Myr) 




7.9 


5.1 


4.6 


Abs. factor 




1.621 


2.082 


2.422 


log Q(H) (s-i) 


Obser. 


52.59 


52.20 


52.34 




Model 


52.59 


52.26 


52.29 


EW(H/3) (A) 


Obser." 


132 


153 


167 




Model 


127 


158 


164 


log Filling factor 




-2.52 


-2.02 


-2.22 


logU 




-2.84 


-2.54 


-2.62 


log Dust-to-gas ratio 




-1.96 


-1.97 


-1.81 


Av 




0.32 


0.49 


0.55 



''Corrected for the underlying stellar population 



modelled values is excellent, both for the number of ionizing 
photons and the EW(H/3). 

In Fig. 5 we show a comparison between the four mea- 
sured electron temperatures in each of the three knots and 
the values predicted by the models. As we can see, the best 
agreement is found for Te([Olll]). In Te([Slll]), a good agree- 
ment is found only for knot A. The model temperatures are 
higher for Te([Sn]) and lower for Te([On]) than the derived 
from the measured line intensities. In Fig. 6, we see the same 
comparison for the four respective ionic abundances. As in 
the case of electron temperatures, the agreement between 
O^^ abundances derived from the observations and found 
by the models is excellent, while in the case of 0+ and S^+, 
only deviations not larger than 0.1 dex are found. The most 
evident deviation is found for the values of S'''/H"'" which 
are higher in the direct measurements than in the model by 
0.3 dex in average for the three knots. 

We have calculated the total abundances of all the mea- 
sured ions, now using these models, taking into account 
when required, the unseen ionization stages of each element, 
using the appropriate model predicted ICF for each species. 
The predicted ICFs for O, S, N, Ar and Ne and the to- 
tal abundances obtained are listed in Table 6. A discussion 
about the differences between the ICFs calculated by these 
models and those obtained from the most commonly used 
formulae is found in Appendix A of Perez-Montero et al. 
(2010). 




a 1.2 




1.0 1.2 1.4 

T ([OIII]) (derived) 



1.0 1.2 1.4 1.6 

T ([Oil]) (derived) 




1.2 1.4 1.6 

T ([SIII]) (derived) 



0.8 1.0 1.2 

T ([SII]) (derived) 



Figure 5. Measured electron temperatures for the three knots, 
vs. values predicted by the photoionization models described in 
the text. The symbols are: grey squares, knot A; red circles, knot 
B, and yellow triangles, knot C. Temperatures are in units of 
10* K. 



4 DISCUSSION 

4.1 Gaseous physical conditions and element 
abundances 

4.1.1 Densities and temperatures 

Four electron temperatures - Te([Olll]), Te([Oll]), Te([Sni]) 
and Te([Sn])- have been estimated for the star-forming 
knots of SDSS J1657. The good quality of the data allows us 
to reach high precision, with rms errors of the order of 2%, 
5%, 6% and 6% in knot A for Te([Oni]), T<,([On]), Te([Sin]), 
and Te([Sll]), respectively. For the faintest knots, B and C, 
the fractional errors are slightly higher, 3%, 8%, 6% and 
17%, respectively. 

The star-forming regions show temperatures within a 
relatively narrow range, between 12000 and 14800 K for 
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Figure 6. Modelled vs. derived oxygen and sulphur ionic abun- 
dances (see description in the text). Knot A: grey squares; knot 
B: red circles; knot C: yellow triangles. 

Te([Olll]). It is worth remembering that the adopted se- 
lection criteria for SDSS J1657 was high H/3 flux and large 
equivalent width of Ha, which tend to render objects with 
abundances and electron temperatures close to the median 
values shown by Hll galaxies. Although these criteria ap- 
plied to the main knot (the SDSS spectrum), we find similar 
electron temperatures for all the regions. To our knowledge, 
there is no previously reported Te([Olll]) for this galaxy in 
the literature. The estimated [Olll] temperature for knot B 
is very similar to that of knot C; both are higher than the 
Te([Oin]) for knot A by about 2000K. At same the time, 
although differences in Te([Sni]) among the three knots are 
much larger, being this temperature 1900 K larger in knot B 
than in knot A, and 1600 K lower in knot C than in knot A, 
these deviations are still compatible within the errors with 
the empirical relation found in Paper I between Te([Oni]) 
and Te([Sni]) for a heterogeneous sample of Giant HII Re- 
gions and Hll galaxies. 

4-1.2 Chemical abundances 

The abundances derived for the three knots using the direct 
method show the characteristic low values found in strong 
line Hll galaxies (Terlevich et al. 1991; Hoyos & Di'az 2006). 
These values are in the range of 12-|-log(0/H) =7.78-7.99, 
in very good agreement with what is found from the pho- 
toionization models discussed above, ranging between 7.80- 
8.02. The data presented in this paper is of high quality and 
the mean error values for the oxygen and neon abundances 
are 0.05 dex, 0.12 for sulphur and 0.06 for argon. Knots B 
and C show a similar value of 12-|-log(0/H), while Knot 
A is almost 0.2 dex higher. This difference is greater than 
the estimated observational errors, and is similar (or even 
smaller) to what is found in other works with spatial reso- 
lution of knots that belong to Hll galaxies or Blue Compact 
Dwarf (BCD) galaxies (see e.g. Izotov et al. 1997; Vflchez & 
Igiesias-Paramo 1998; Kehrig et al. 2008; Cairos et al. 2009; 



Perez- Montero et al. 2009; Garci'a-Benito et al. 2010). How- 
ever, in general, these differences were attributed to the ob- 
servational uncertainties (pointing errors, seeing variations, 
etc.) or errors associated to the reddening correction and flux 
calibration, and the oxygen abundance variations were not 
assumed as statistically significant, concluding that there 
is a possible common chemical evolution scenario in all of 
them. There are even greater differences when comparing 
the estimated abundances of the individual knots with those 
derived from the integrated spectra of the galaxies. For in- 
stance, Cairos et al. (2009) found for the integrated spec- 
trum of Mrk 1418 a lower value of direct oxygen abundance 
by about 0.35 dex (equivalent to a factor of 2.2) than for 
knots 1 and 2 of that galaxy. They pointed out that while 
this variation could reflect a real abundance difference in 
different scales (kpc-sized aperture for the integrated spec- 
trum and sizes of the order of 100 pc for individual Hll re- 
gions), it may also be due to relatively large measurement 
uncertainties for the weak [Olll] auroral emission line. For- 
tunately, our data are not affected by pointing errors and 
seeing variations, and the other observational uncertainties 
have a second order effect, since TWIN is a double beam 
longslit spectrograph that simultaneously acquire all the ob- 
served spectral range. Likewise, the errors associated with 
the measurements of the weak auroral emission lines are rel- 
atively small, specially for [Olll]. 

The logarithmic N/0 ratios found for SDSS J1657 us- 
ing the direct method are -1.23 ±0.11, -1.35 ±0.16 and - 
1.36 ±0.15 for knots A, B, and C, respectively. The derived 
values are on the high log(N/0) side of the distribution for 
this kind of objects (see left-hand panel of Fig. 6 of Pa- 
per II). The logarithmic values of this ratio found for the 
three knots using photoionization models are slightly lower, 
-1.34 ±0.07, -1.47 ±0.11 and -1.43 ±0.11, respectively, al- 
though similar within the errors. However, the derived val- 
ues present a larger uncertainty than the values of total oxy- 
gen abundance, so a deflnite conclusion can not be extracted 
about the homogeneity of this ratio. Anyway, it is quite sug- 
gestive to find a larger N/0 ratio in the brightest knot, which 
has the larger metallicity. The N/0 ratio is directly related 
to the chemical history of galaxies, as these two elements 
have different nucleosynthetic origin, so this difference can 
support to some extent the idea of a different chemical evo- 
lution in the three knots of this galaxy. 

The log(S/0) ratios found are quite similar for Knots 
A and B (-1.53 and -1.57 respectively) and higher for Knot 
C (-1.32), but all three are almost equal within observa- 
tional errors (~ 0.12 in average). The values derived using 
the photoionization models follow the same trend as those 
estimated using the direct method (-1.65, -1.69, and -1.42 re- 
spectively). They are all slightly lower than the solar value, 
log(S/O)0 =-1.36 (Grevesse & Sauval 1998). On the other 
hand, the logarithmic Ne/O ratio is remarkably constant, 
with a mean value of 0.75 (0.72 from the photoionization 
models), despite the differences in oxygen abundance be- 
tween knot A and knots B and C. They are consistent with 
solar, log(Ne/0) =-0.61^. The Ar/O ratios found (which are 
almost the same using the direct method and the photoion- 



* Oxygen from Allcnde-Pricto et al. (2001) and neon from 
Grevesse & Sauval (1998). 
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Figure 7. Left panel: log(0+/0^+) vs. log{S+/S'^"'") using the direct method and the photoionization models (filled yellow diamonds 
and red triangles), the objects studied in Paper I and Paper II (blue circles) and Hll galaxies from the literature (open squares). Diagonals 
in this diagram correspond to constant values of 77. Right panel: log([On]/[Oin]) vs. log([Sll]/[Slll]), symbols as in left panel. Diagonals 
in this diagram correspond to constant values of rj'. 



ization models) show a very similar value for Knot A and 
B, while Knot C has a ratio higher by 0.2 dex. The mean 
value is consistent with solar, log(Ar/0) =-2.29 ', within the 
observational errors. 

Finally, the derived helium abundances are the same for 
the three knots within observational errors. 



4.2 Ionization structure 

The ionization structure of a nebula depends essentially on 
the shape of the ionizing continuum and the nebular ge- 
ometry and can be traced by the ratio of successive stages 
of ionization of the different elements. With our data it is 
possible to use the O'^ /O^'^ and the S^/S^"*" to probe the 
nebular ionization structure. In fact, Vilchez & Pagel (1988) 
showed that the quotient of these two quantities that they 
called "softness parameter" and denoted by rj is intrinsi- 
cally related to the shape of the ionizing continuum and 
depends on geometry only slightly. An insight into the ion- 
ization structure of the observed objects can be gained by 
means of the 77 diagram (see Paper I). 

In Fig. 7, left panel, we show the relation between 
log(0+/0^+) and log(S+/S^+) derived using the direct 
method and the photoionization models for the knots of 
SDSS J1657 (filled yellow diamonds and red triangles, re- 
spectively), the objects studied in Paper I and Paper II (blue 
circles) and Hll galaxies (open squares) from the literature 
(see description and references in Paper II). In this diagram 
diagonal lines correspond to constant values of the rj pa- 
rameter which can be taken as an indicator of the ionizing 
temperature (Vilchez & Pagel 1988). Hll galaxies occupy the 
region with log rj between -0.35 and 0.2, which corresponds 
to high values of the ionizing temperature. As noticed in 
Paper II, Knot A shows a very low logarithmic value of 



^ Oxygen from AUcnde-Pricto et al. (2001) and argon from 
Grevesse & Sauval (1998). 



rj, -0.6. Knots B and C present very similar values. These 
objects, however, have the [On] AA 7319,25 A lines affected 
by a set of atmospheric absorption lines. We treated the 
problem (as is custom) by dividing the spectra by a tel- 
luric standard. Unfortunately no previous data of this ob- 
ject exists apart from the SDSS spectrum of Knot A. While 
the agreement between the [Olll] lines measured in both the 
SDSS spectrum and ours is good, the [On] lines measured on 
the SDSS spectrum provide a Te([On]) = 1.23 ± 0.21, lower 
than the one derived from our data. This lower tempera- 
ture would increase the value of 0+/0^+ moving the cor- 
responding data point upwards in the left panel of Fig. 7. 
This would be consistent with the position of the object 
in the right panel of the figure showing log([On]/[Oni]) vs. 
log([Sll]/[Slll]), which does not require explicit knowledge 
of the line temperatures involved in the derivation of the 
ionic ratios, and therefore does not depend on the method 
to derive or estimate these temperatures. The right panel 
in Fig. 7 shows the purely observational counterpart of the 
left panel. In this diagram diagonal lines represent constant 
values of log ry ' = log [( [Oil] / [Olll] ) / ( [Sll] / [Sin] )] . 

Another possible explanation for the differences be- 
tween these two diagrams can be obtained by inspecting 
the position of the models described above in relation to the 
observational points (red triangles in the figures). Models 
with a thick shell geometry and a constant density predict 
higher Te([Sll]) and, hence, lower S"'"/H"'', which causes the 
rj parameter to be higher than the values estimated from 
the measurements. In fact, the difference between the ionic 
abundances derived by the direct method and those pre- 
dicted by the models can be also seen in Fig. 6. On aver- 
age there is a difference of 0.23 dex for O^^/H^, 0.07 dex for 
02+/H+, 0.67 dex for S+/H+ and 0.13 dex for S'^+/R+. This 
effect, already pointed out by Perez-Montero et al. (2010), 
could be a consequence of an outer shell of cold diffuse ion- 
ization structure in these objects with an extra emission of 
[Sll] which contributes to their integrated spectrum. The 
agreement is better for the rj' diagram, being knot B the 
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Figure 8. The oxygen abundances and their uncertainties for each observed knot of SDSS J1657 (blue solid Une rectangles), as derived 
using different empirical calibrators. From left to right top panels: O23, N2, and 03N2 and bottom ones: S23, S23/O23, and Ar303. The 
dash (black) and dotted (red) line rectangles represent the abundances and their uncertainties as derived from the direct method and 
the photoionization models, respectively. 



most discrepant, which is consistent with the differences be- 
tween the observed and modelled lines (see Fig. 4), where 
the higher difference corresponds to the [Sll] lines. 

In both diagrams, 77 and 77', the three knots of SDSS 
J1657 present a very similar ionization structure, showing 
almost the same values within the observational errors. This 
implies that the equivalent effective temperatures of the ion- 
ization radiation field are very similar for all the knots, al- 
though we find some small differences in the ionization state 
of the different elements. 



4.3 Chemical abundances from empirical 
calibrators 

The emission line spectra of the three star-forming knots in 
J1657 are very similar, implying similar values for ioniza- 
tion parameter, ionization temperature, and chemical abun- 
dances. We derived the ionization parameters from the [Oil] 
to [Olll] lines ratio according to the expression given in Di'az 
et al. (2000) . The logarithmic ratio is similar in all the knots 
ranging from -2.47 for knot B to -2.65 for knot A. 

The different strong-line empirical methods for abun- 
dance derivations, which have been widely studied in the 
literature, are based on directly calibrating the relative in- 
tensity of some bright emission lines against the abundance 
of some relevant ions present in the nebula (see e.g. Garcfa- 
Lorenzo et al. 2008; Cairos et al. 2009; Garci'a-Benito 2009; 
Garcia- Benito et al. 2010, and references therein). For the 
case of oxygen, we take the calibrations studied by Perez- 
Montero & Dfaz (2005), who obtain different uncertainties 
for each parameter in a sample of ionized gaseous nebulae 



with accurate determinations of chemical abundances in the 
whole range of metallicity. 

In Fig. 8, we show the total abundances as derived 
from several strong-line empirical methods (with their cor- 
responding errors estimated taking into account the errors 
of the line intensities and also the errors given by the cali- 
brations of the empirical parameters) and the oxygen abun- 
dances calculated from the electron temperatures measured 
using the direct method and those estimated from the pho- 
toionization models for each knot. 

Among the available strong-line empirical parameters 
we studied the O23 parameter (also known as R23 and orig- 
inally defined by Pagel et al. (1979) and based on [On] and 
[Olll] strong emission lines). This parameter is character- 
ized by its double-valued relation with metallicity, with a 
very large dispersion in the turnover region. According to 
the values measured, we used the McGaugh (1991) calibra- 
tion for the lower branch. For knots B and C, this calibrator 
fails to predict the value obtained with the direct method, 
overestimating the oxygen abundance, although the derived 
values are very similar if we take into account the obser- 
vational errors and the large spread in the empirical O23 
diagram (see Fig. 3 of Perez-Montero & Dfaz 2005). 

The N2 parameter (defined by Storchi-Bergmann et al. 
1994) is based on the strong emission lines of [Nil]. It remains 
single-valued up to high metallicities in its relation to oxy- 
gen abundance, and it is almost independent of reddening 
and flux calibrations. Nevertheless, it has the high disper- 
sion associated to the functional parameters of the nebula 
(ionization parameter and ionizing radiation temperature) 
and to N/O variations. We used the empirical calibration 
of this parameter from Denicolo et al. (2002) to derive the 
oxygen abundance in the three star-forming knots of this 
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galcixy. We can see in Fig. 8 that N2 behaves similarly to 
O23 in predicting the abundances. 

The parameter 03N2, defined by AUoin et al. (1979), 
depends on strong emission lines of [Olll] and [Nil]. We used 
the calibration due to Pettini & Pagel (2004) and, as we can 
see in Fig. 8 it has a very similar behaviour to that of N2. 

The S23 parameter was defined by Vilchez & Esteban 
(1996) and is based on the strong emission lines of [Sll] 
and [Slll] . The cahbration by Perez-Montero & Dfaz (2005) 
yields comparable oxygen abundances for the three observed 
knots that are in turn in very good agreement with the abun- 
dances derived using the direct method for knots A and B, 
and slightly higher for knot C, but still consistent within the 
errors. 

The ratio of the S23 and O23 parameters as S23/O23 
(Dfaz & Perez-Montero 2000) is a parameter that increases 
monotonically with the oxygen abundance up to the overso- 
lar regime and is very useful to study variations over wide 
ranges of metallicity {e.g. disks). We applied the calibration 
from Perez-Montero & Diaz (2005) and found a good con- 
cordance with the values determined by the direct method. 

The ArsOs parameter, defined and calibrated by 
Stasihska (2006) as the ratio of [Arm] A 7136 A to 
[Olll] A 5007 A lines, predicts slightly higher values than the 
S23 parameter and in better agreement with those derived 
directly. 

Different strong line empirical metallicity calibrators 
are commonly used to estimate the oxygen abundances in 
objects for which direct derivation of electron temperatures 
is not possible. However, as illustrated in Fig. 8, different 
empirical calibrations give results that are not in complete 
agreement with direct measurements, and the goodness of 
the result even changes between knots when using the same 
calibrator. Perez-Montero et al. (2009) found a very simi- 
lar behaviour for the star- forming knots of IIZw71, except 
for the Ar303 parameter, which is a better estimator of the 
oxygen abundances for the knots in SDSS J 1657. For two 
knots in Mrk 1418, Cairos et al. (2009) derived the oxy- 
gen abundances from the observed [Olll] A 4363 A auroral 
emission line, and compared them with those abundance 
values derived using the N2 and 03N2 empirical parame- 
ters, and found a similar discrepancy, with the empirically 
derived values being slightly larger. Clearly, more observa- 
tions and direct abundance estimations are needed in order 
to improve calibrations and truly understand the origin of 
the observed dispersions and discrepancies (see discussion 
in Paper I about the dispersion and precision of the S23 
parameter) . 

4.4 The stellar population 

The study of the stellar content of our objects was car- 
ried out using the STARLIGHT code, which calculates the 
combination of stellar libraries and the extinction law that 
reproduces the spectral energy distribution, to derive the 
properties of the stellar population in each of the knots. 
STARLIGHT fits an observed continuum spectral energy 
distribution using a combination of multiple simple stellar 
populations (SSPs; also known as instantaneous bursts) syn- 
thetic spectra using a minimization procedure. For con- 
sistency with the photoionization models used to model the 
gas of knot A in Perez-Montero et al. (2010), we have used in 



Table 8. Values of the extinction, total stellar mass, and fraction 
of the mass in stars younger than 10 Myr for each knot of SDSS 
J1657 in the best-fit model using STARLIGHT. 



ID 


A(V) 


c(H/3) 


M, 




f(age < 10 Myr) 




(mag) 








(%) 


Knot A 


0.18 


0.08 


280 


1.6 


0.58 


Knot B 


0.00 


0.00 


100 


0.2 


0.21 


Knot C 

M^4-^ TV /T 


0.00 


0.00 


140 


0.4 


0.28 



log(Age) 

Figure 9. Histogram of the distribution in visual light of the 
most probable stellar population models fitted by STARLIGHT 
for knots A, B and C, as labelled. 



the fitting a synthetic stellar population obtained using Star- 
burst99 (Leitherer et al. 1999; Vazquez & Leitherer 2005) 
with the Geneva stellar evolutionary tracks for continuous 
star formation with high mass loss (Meynet et al. 1994), the 
Kroupa Initial Mass Function (IMF; Kroupa 2002) in two 
intervals (0.1-0.5 and 0.5-100 M©) with different exponents 
(1.3 and 2.3, respectively), the theoretical wind model (Lei- 
therer et al. 1992), with the model atmospheres from Smith 
et al. (2002), and the stellar cluster metallicity being the 
closest to the nebular one, Z = 0.004 (=1/5Z0, see Paper 
II and Perez-Montero et al. 2010). The STARLIGHT code 
solves simultaneously the ages and relative contributions of 
the different SSPs and the average reddening. The reddening 
law from Cardelli et al. (1989) is used. Prior to the fitting 
procedure, the spectra were shifted to the rest-frame, and 
re-sampled to a wavelength interval of 1 A in the entire wave- 
length range between 3500 A and 9000 A by interpolation 
conserving flux, as required by the program. Bad pixels and 
emission lines were excluded from the final fits. 

In Fig. 9 we show the age distribution of the visual light 
fraction for the individual knots. All of them present a very 
young stellar population with ages around 10 Myr, responsi- 
ble for the ionization of the surrounding gas. Practically all 
the mass of the knots seems to come from a very old popula- 
tion of about 8.3 Gyr. However, this result is puzzling given 
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that no absorption metal lines characteristic of old stellar 
populations, such as MgiiA5173A and/or CaiiH A 3933 A, 
Can K A 3968 A, and Can Triplet AA8494, 8542, 8662 A, are 
detected in the spectra. The estimated total stellar mass, 
the mass of the stellar population with an age younger than 
10 million years, responsible for the ionization of the gas, 
and the fraction of this last one with respect to the former, 
are listed in Table 8 for the three knots, together with the 
internal extinction, A(V) (as given by STARLIGHT) and 
the reddening constant [c(H/3)] estimated by the model. No 
aperture correction factors have been taken into account for 
the Ha luminosities, due to the compact nature of the ob- 
jects. Indeed, the discrepancy factors between our Ha lumi- 
nosities for Knot A and those measured in the SDSS catalog 
using a 3 arcsec fiber is no larger than 1.3. 

We do not expect to find the same values of extinction 
in the gas and the stellar population. In this galaxy, although 
the extinction is larger in knot A than in B and C (the op- 
posite from the derived values using the Balmer decrement) 
the general result is consistent with low extinction. 

The mean ages and stellar mass fraction (with respect 
to the total stellar mass of the corresponding cluster) of the 
ionizing stellar population fitting by STARLIGHT for each 
knot is very similar, in very good agreement with the ion- 
ization structure derived in §4.2. Since the age distribution 
of the ionizing population of the different knots seems to 
be similar, this could indicate a common evolutionary stage 
of this population which is probably related to a process of 
interaction with a companion galaxy that triggered the star 
formation in the different knots almost at the same time. 
In Fig. 1 we can appreciate a tail of diffuse emission gas 
observed to the south of knot A, with a few weak and blue 
knots probably hosting star formation. To the north-east 
there is a bright and redder galaxy. We need more spec- 
troscopic data to disentangle whether these different stellar 
systems are related or it is only a projection effect. The rel- 
atively small differences in the derived metallicities for the 
star-forming knots of SDSS J1657 could indicate a slightly 
different chemical history of these knots. 



4.5 Ionizing stellar populations 

Some properties of the emission knots can be obtained from 
the measured Ha flux, such as Ha luminosity, number of 
ionizing photons, mass of ionizing stars and mass of ionized 
hydrogen (see e.g. Diaz et al. 2000). The observed Ha flux 
was corrected in each knot for reddening using the values 
of the reddening constants, c(H/3), given in Table 3. We as- 
sume a distance to SDSS J1657 of 161.2 Mpc (Mould et al. 
2000). Besides, we have used the dust absorption factors, fd 
derived using photoionization models to correct these quan- 
tities. These factors are independent of reddening correction 
because they affect above all the number of ionizing photons 
emitted by the stellar cluster. As we can see in Table 7, these 
factors are larger in knots B and C than in knot A, consis- 
tent with the differences found in the reddening as derived 
by measuring the Balmer decrement. These differences in 
reddening are not surprising attending to the different ages 
of the ionizing clusters found in the same photoionization 
models (see §3.4), which are respectively 7.9, 5.1 and 4.6. Al- 
though the metallicity in knots B and C is lower, as the star 
formation started later in these knots, the ionizing photons 



Table 9. Derived properties of the observed knots using the ex- 
tinction and dust absorption corrected measured Ha fluxes. 



ID 


F(Ha) 


L(Ha) 


Q(Ho) 




M(Hn) 


SFR 


Knot A 


28.4 


8.3 


6.3 


3.6 


2.1 


0.697 


Knot B 


15.2 


4.8 


3.3 


1.5 


3.5 


0.375 


Knot C 


23.0 


7.0 


5.3 

-1 


2.2 

-■2 


5.3 


0.564 



10**^ ergs ^, ionizing photons in 10^^ photons s ^, masses 
in 10'^ Mq, and SFR in Mq yr'^ 



have not already broken the dust cocoon usually shroud- 
ing these very young star formation knots. The derived and 
corrected values are given in Table 9. 

The derived values of the masses of the ionizing clusters 
(i.e. with an age younger than 10 Myr) can be compared 
with those provided by the STARLIGHT fit. For knot A 
STARLIGHT gives a slightly lower value than the one we 
derive by a factor of 2, while for knot B and C this factor 
is larger, about 7 and 5, respectively. If we take the masses 
derived from the Ha fiuxes, which are, in principle, lower 
limits since a possible escape of photons is not taken into 
account, and we use the calculated proportions of young to 
total mass given in Table 8 we obtain total masses of about 
7 X 10** M0 for each knot. 

The weight of the underlying stellar population in the 
visual light is larger in knot A than in knots B and C. The 
corrected EW(H/3) is 12% larger in knot A but less than 1% 
in the other two knots. 

The star formation rate (SFR) for each knot was de- 
rived from the Ha luminosity using the expression given by 
Kennicutt (1998), SFR = 7.9 x 10"*^ x L(Ha). The derived 
values are also given in Table 9, taking into account dust 
absorption factors obtained from our models as discussed in 
previous sections. The SFR obtained range from 0.375 Mq 
yr"^ for knot B to 0.697 Mq yr"^ for knot A. 



4.6 Kinematics of the knots 

In Fig. 10 we show the longslit spectrum between 5030 
and 5220 A (upper panel) and between 6690 and 7110 A 
(middle panel). These spectral ranges contain H/3 and 
[Qui] 4959,5007 A, and Ha, [Nil] 6548,6584 A, Hei6678A 
and [Sii] 6717,6731 A emission lines, respectively. All these 
lines show a similar behaviour in the position- velocity (wave- 
length) plane, with the typical helical shape of the rotation 
curve expected for a rotating disc. Knot B is redshifted while 
Knot C is blueshifted with respect to the systemic velocity 
of knot A. We have analyzed the differential radial velocity 
field along the slit using the Ha and [Olll] emission lines. 
The radial velocities derived for each pixel (0.56 arcsecs) are 
displayed in the lower panel of Fig. 10. We have also plot- 
ted the Ha spatial profile. The velocities derived from both 
emission lines are in good agreement, giving the same value 
within the observational errors, which are about 20km/s, 
taking into account the errors in the wavelength calibration 
and a single Gaussian fitting. Larger differences (~ 30km/s) 
between the derived velocities using the different emission 
lines have been found around the secondary knots, B and 
C, but they are also within the observational errors. From 
the data for Knot A, we have estimated a redshift of 0.0383 
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Figure 10. Longslit spectrum between 5030 and 5220 A (upper 
panel) and between 6690 and 7110 A (middle panel); the main 
emission lines and the knots are labelled. The derived rotation 
curve of SDSS J1657 is shown in the lower panel. Green circles 
represent the radial velocities derived from Ho and yellow dia- 
monds from [Oin]5007A for each pixel (0.56 arcsec) along the 
slit. Red squares and blue upward triangles represent the radial 
velocities derived from Ha and [Olll] using the extracted spectrum 
for each knot. The dashed line represents the spatial distribution 
of the Ha flux to be compared with the rotation curve. 



(vr 11500 km/s), in very good agreement with the value 
given by the SDSS (see Table 1). 

The rotation curve is basically linear, with some per- 
turbations. Also the weak pure emission line knot located 
between Knot A and C seems to follow a circular rotation 
curve. There is a depression of about 50km/s in the radial 
velocity curve located between this pure emission knot and 
Knot C. In principle, we could attribute this deviation from 
the circular motion to a poor signal from the emission lines 
(see the spatial profile and the longslit spectra. Figs. 2 and 



10, respectively). However, it is not the position with the 
lowest emission, and the other zone around the minimum 
emission does not show a similar deviation of the radial ve- 
locity field from the circular motion. Moreover, the radial 
velocities in this depression derived from the Ha and [Oin] 
lines are in good agreement showing a similar behaviour, 
which leads us to suppose that this depression is a real de- 
viation from the circular motion. This could be due to the 
presence of an expanding bubble or shell (or superposition of 
more than one) of ionized gas approaching us with respect to 
the systemic velocity of the galaxy defined by the velocity of 
the brightest and most massive star-forming cluster. Knot 
A. In this case, we could only see the expanding material 
which is moving in our direction, since we detect a system- 
atic variation of the radial velocity, but not a broadening 
of the emission line with a systemic velocity in accordance 
with that expected for the circular motion. Another pos- 
sibility would be that the ionized gas related to this weak 
emission area is affected by a possible interaction with a tail 
of emission gas located at the south of the Knots or with the 
redder galaxy observed to the north-east (see Fig. 1 and dis- 
cussion above). This interaction could be responsible for the 
deviation from circular motion, as well as the fact that the 
ionizing population of the different star-forming knots seem 
to have similar ages (see §4.4). It has been noticed by Cairos 
et al. (2009) that several works show a substantial fraction 
of dwarf galaxies that show recent episodes of star formation 
have optical faint and gas-rich low-mass companion galax- 
ies (e.g. Taylor et al. 1993, 1995, 1996; Noeske et al. 2001; 
Pustilnik et al. 2001). Even in the case of apparent isolated 
galaxies there are clues that favour the interpretation that 
interactions play a substantial role in the evolution of these 
individual systems (Johnson et al. 2004; Ostlin et al. 2004; 
Bravo- Alfaro et al. 2004, 2006; Gumming et al. 2008; Cairos 
et al. 2009). 

Assuming that we see the galaxy edge-on, with a radial 
velocity Vr = 100 km/s given by the difference between Knot 
A and the furthest point where we can measure the emission 
lines with good enough S/N, and considering an optical ra- 
dius of 14 arcsec for this point (equivalent to Rgpt = 11 Kpc 
at the adopted distance of 161.2 Mpc for SDSS J1657, with 
a scale of 782 pc arcsec"^) we have estimated a dynamical 
mass inside this radius of 2.5 x 10^° Mq. 



Since the three star-forming knots clearly dominate the 
luminosity of the galaxy (see Fig. 1), we have derived their 
blue magnitudes from the SDSS values (see Table 1). We 
adopt these values as representative of the total light, al- 
though we know that they are (together) an upper limit to 
the total magnitude of the galaxy since we are not taking 
into account the light outside the 3 arcsec SDSS fiber (see 
Fig. 1). We have obtained the B magnitude from the Sloan 
g and r photometric magnitudes using the transformation 
equation from Chonis & GaskeU (2008), B = g-f 0.327x (g- 
r)-f0.216, given B^ = 17.15, Bb=20.75, and Bc = 19.67. We 
have obtained a total blue luminosity Lg = 6.2x 10® L0, 
which gives an upper estimation of 4 for M/Lg , in agreement 
to what is found by Faber & Gallagher (1979) for irregular 
galaxies. 
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5 SUMMARY AND CONCLUSIONS 

The star formation processes in Hll galaxies are known to oc- 
cur in low density environments. Spectrophotometric obser- 
vations of SDSS J1657 were carried out to study the physical 
properties of its individual bursts of star formation. In the 
three individualized star-forming regions the electron den- 
sities were found to be well below the critical density for 
coUisional de-excitation. 

We extracted information on the three knots of SDSS 
J1657, labelled A, B and C. For all of them we measured 
four line temperatures: Te([Olll]), Te([Sni]), Te([Oll]), and 
Te([Sll]), reaching high precision, with rms fractional errors 
of the order of 2%, 5%, 6% and 6%, respectively, for Knot A, 
and slightly worse 3%, 8%, 6% and 17%, respectively, for the 
fainter knots B and C. The [Olll] temperature of knot A was 
found to be about 2000 K lower than for the other two knots. 
For the [On] temperature the estimated values for the three 
knots are the same within the errors. The [Slll] temperatures 
show a difference of 3500 K between the values derived for 
Knots B and C, and the estimated value for Knot A falls 
in between. However, these large differences are compatible 
with the empirical relation between Te([Oni]) and Te([Sni]) 
found in Paper I, within the observational errors and the 
empirical dispersion of that relation. Within observational 
errors, the [Sll] temperature derived is similar for the three 
knots. 

The temperature measurements allowed the direct 
derivation of ionic abundances of oxygen, sulphur, nitrogen, 
neon and argon. The total abundances of these species are in 
the same range of metallicities measured in Hll galaxies, be- 
tween 7.78 and 7.99, with estimated errors of about 0.05 dex. 
Knots B and C show similar abundances, while the value for 
Knot A is about 0.2 dex higher. This behaviour is mirrored 
by the N/O ratio. Knots A and B have similar Ar/O while 
Knot C is the one with a different value. Within obseva- 
tional errors, S/0 is almost the same for the three knots. 
The Ne/0 ratio is remarkably similar for all three regions. 

We have studied the underlying and ionizing stellar pop- 
ulations for the three regions and modelled the properties 
of the emitting gas using a photoionization code. The es- 
timated electron temperatures and ionic abundances using 
the direct method are well reproduced by the models except 
for those of S^/H+ . This could be due to the presence of 
diffuse gas in these star-forming regions, which is not taken 
into account in the models. Regarding the hardness of the 
radiation field, model predictions agree with observations 
when the softness parameter rj which parametrizes the tem- 
perature of the ionizing radiation, is expressed in terms of 
emission lines intensity ratios. However, there is a disagree- 
ment when ionic abundance ratios are used instead. This is 
probably caused by the overestimate of the electron temper- 
ature of S"*" by the models and the corresponding underes- 
timate of the S^/S'^"'' abundance ratio rather than by the 
existence of an additional heating source in these knots. 

We have also estimated the total oxygen abundances by 
means of different strong-line empirical parameters. In all 
cases, the estimated abundances are consistent with those 
derived by the direct method with the parameter involving 
the sulphur lines providing the closest fit. The rest of the 
parameters slightly overestimate the oxygen abundance. 

The star formation history for the three knots derived 



from the fitting of multiple simple stellar populations using 
STARLIGHT are remarkably similar, with an old popula- 
tion of about 8 Gyr presenting more than the 99% contri- 
bution to the mass fraction. This result is somewhat unex- 
pected and data on its surface brightness distribution would 
be very valuable to explore further this finding. Since the age 
distributions of the ionizing population among the different 
knots of SDSS J1657 seem to be similar, this could indicate 
a common evolutionary stage of this population which is 
probably related to a process of interaction with a compan- 
ion galaxy that triggered the star formation in the differ- 
ent knots at about the same time. Besides, this interaction 
could be responsable for the deviation from the circular mo- 
tion shown by the rotation curve, which is basically linear 
in the region where we can measure the emission lines with 
a good enough S/N ratio. However, this deviation could be 
also related to an expanding bubble or shell of ionized gas 
approaching us with a velocity of 50km/s with respect to 
the predicted velocity from the rotation curve. 

The ionization structure mapped through the use of the 
ri and rj' diagrams derived from our observations shows very 
similar values within the errors for all the knots. This fact 
implies that the equivalent effective temperatures of the ion- 
ization radiation fields are very similar for all the studied 
regions, in spite of some small differences in the ionization 
state of different elements. 

Finally, we have derived a dynamical mass of 
2.5 X W^^ Mq from the rotation curve with a mean value 
of lOOkm/s at an radius of llKpc. The total mass of the 
young clusters derived for the three star-forming knots us- 
ing the Hq luminosities is 7.3 x 10^ Mq, making up a small 
fraction of the total dynamical mass, about 0.03 %. We have 
estimated an upper limit of about 4 for the ratio M/Lg. 
Data on the surface brightness distribution using broad and 
narrow band images would be very valuable to explore fur- 
ther this result. 
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